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We show analytic results for the irreducible charge and spin susceptibilities, xo('^,Q), where Q 
is the momentum transfer between the nodes in d-wave superconductors. Using the BCS theory and 
a circular Fermi surface, we find that the singular behavior of the irreducible charge susceptibility 
leads to the dynamic incommensurate charge collective modes. The peaks in the charge structure 
factor occur at a set of wave vectors which form an ellipse around Q-n- = (tt, tt) and Qo = (0, 0) 
in momentum space with momentum dependent spectral weight. It is also found that, due to 
the non-singular irreducible spin susceptibility, an extremely strong interaction via random phase 
approximation is required to support the magnetic peaks near Q^. Under certain conditions, the 
peaks in the magnetic structure factor occur near Q = (tt, 7r(l ± 5)) and (7r(l ± 5), tt). 

PACS numbers; 74.20.-z, 74.25.Nf, 78.70.Nx 

The neutron scattering measurements provide the direct information about the wavevector and frequency depen- 
dence of the dynamic spin susceptibihty. On the other hand the inelastic x-ray or electron scattering can measure the 
dynamic charge susceptibility. These informations are particularly important in the cuprate superconductors because 
of the intimate interplay between the spin and charge dynamics which may be related to the mechanism of the high 
temperature superconductivity. 

The neutron scattering experiments on La2-xSrxCu04{LSCO) j^J|] found the inelastic incommensurate peaks 
near Q^r = (tTjTt) in the superconducting state as well as in the normal state The collective modes occur at 
Q = (7r,7r(l ± S)) and (7r(l ± S),n) for low frequencies. Recently, a lot of effort was put forward to investigate the 
spin dynamics in Y Ba2Cu^Oj^y{Y BCO) In earlier works, the incommensurate peaks were reported in the scan 

along the zone diagonal direction |^ . Later it was found that the locations of the higher intensity peaks are the same 
as those of LSCO; (tt, 7r(l ± 5)) and (7r(l ± (5), tt) [|. 

On the theoretical front, there have been two different approaches to explain the incommensurate peaks. The nu- 
merical studies of one- and three-band Hubbard models found that the incommensurate peaks occur as the dynamical 
response of a spatially homogeneous system with a nearly nested Fermi surface. It was also pointed out that the 
nesting peaks are at (7r(l ± 5), tt{1 ± 5)) and (7r(l ± 5), 7r(l =F 5)) in YBCO due to the 45 degree rotation of the Fermi 
surface. |0] From a different viewpoint, it was proposed that the magnetic incommensurate peaks are induced by the 
dynamic charge stripes in a spatially inhomogeneous system. |^ The static charge ordering can occur in the relatively 
bad metal due to the pinning of the charge stripes. On the other hand, the dynamic fluctuations of the charge stripes 
result in the phase coherence of the superconducting state. 

Recently, Tranquada et al found, in the neutron diffraction measurements, the incommensurate elastic magnetic 
peaks at Q = (tt, tt{1±S)) and {tt{1±6), tt) in Lai.Q-xNdoASrxCuO^. They interpreted these peaks as the evidence 
for the static stripe order of charge and spins. It was also shown that the superconducting transition temperature 
increases as the peak splitting parameter, S, increases. This may imply that the charge and spin stripes are intimately 
related to the superconductivity. 

These elastic and inelastic neutron scattering experments motivated us to investigate the incommensurate charge 
and spin fluctuations in the superconducting state. In particular, since the inelastic experiment on YBCO found the 
incommensurate peaks at the same positions as those of LSCO, these may not be related to the details of the Fermi 
surface shape. Thus, as a first step, it is worthwhile to understand the charge and spin dynamics within the BCS 
theory with a simple Fermi surface. This would provide an important information which should be compared with 
any new theoretical proposal-either Fermi surface effect or the charge stripes. 

In this paper, we compute analytically the irreducible charge and spin susceptibilities in the BCS d-wave super- 
conducting state with a circular Fermi surface. We found that the incommensurate charge collective modes occur 
at a set of wave vectors. The wave vectors form an ellipse around Q,r = (tt, tt) in momentum space. The spectral 
weight also depends on the wave vectors. It is found that the incommensurate spin collective mode can occur near 
Q = (tt, 7r(l ± 6)) and (7r(l ± S),'k) only when an extremely strong interaction is assumed via the random phase 
approximation. Using the analytic form of the irreducible spin susceptibility, we also examine the NMR relaxation 
rate 1/Ti for low temperatures and found that (TiT)^^ ^ T'^ / (AEp)^ . We also analyze the quasiparticle lifetime for 
small frequencies at the nodes; 1/r ^ ¥"^0:^ / (AEp)'^ , where V is the interaction. 
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Let us consider the simplest model for the electronic energy with a circular Fermi-surface of radius kp. 

i-2 p 1,2 

a - ek - M = — + — - — , (1) 
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The lowest order charge, Xo^ ^^'^ spin susceptibilities, Xg, for momentum Q and energy w at T = are given by 
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where £'k = + with Ak — A cos 20. The different coherence factors in the charge and spin susceptibilities in 
Eq. (^) come from the fact that the magnetic scattering is odd with respect to the time reversal symmetry while the 
charge scattering is even. 

Let us exmaine the momentum transfer near Q^r = [V^kp, \/2kp) or (tt, tt) shown in Fig. 1. Near the nodes, we 
expand the electronic dispersion relation and the amplitude of the gap. We find the following for |p| + |q| < kp/^/2. 

— -V2vp{px + Py) , Ak ~ -V2A{px - Py)/kp , 

&+Q - \/2vp{px +qx+Py + Qy) , Ak+Q ~ V2A{px +qx~Py- %)/kp . (4) 

Here p = k — kp and q = Q — Qtt, where ko denotes the node of d-wave superconductors. Note also that one obtains 
the same dispersion relations near (tt, tt) in the tight binding model with l/m ~ to? and y^kp = 7r/2a ||lo| , [ll[] , where 
t is the hopping amplitude and a is the lattice spacing. Then, the Dirac dispersion for the quasiparticle is found near 
the nodes; 



Ek = J2EIpI + 2A V- 



Sk+q = ^2El{p+ + q+y + 2A(p_ + , (5) 

where 

P± = {px±Py)/kp , q± = [qx ± qy)/kF . (6) 
Evaluating Eq.^ with Eq.^, the imaginary part of the charge susceptibilities is found for ujq < uj < A, 

where 

UJ^ = 2Elq\ + 2A^ql . (8) 

It can be seen in Eq. (|^) that the charge susceptibility diverges as ljJ approaches ujq and has the following scaling 
form; 



T J V2Epq+ V2Aq^ \ 



where 



Note that the incommensurate collective modes occur with a set of wave vectors. The wave vectors form an ellipse 
around Q,r = (tt, tt) in momentum space for the energy ujq. The essentricity of the ellipse is determined by the 
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ratio of A and Ep. Since Ep is typically larger than A (-^ = 5 ~ 10), the ellipse is elongated along the direction 
perpendicular to the zone diagonal. This result is a consequence of the anisotropic Dirac dispersion relation of the 
quasiparticles near the node. The Eq. (^0|) implies that the spectral weight of the incommensurate collective mode 
depends on the wave vectors. In particular, the spectral weight of the collective mode vanishes at two points in the 
zone diagonal direction. It is found that the charge susceptibility in the random phase approximation, x'^('^j Q) — 
Xo{oj,CI)/[1 + V{lj,CI)xq{i^,Q)], is almost the same as the irreducible susceptibility, Xg(w, Q), for quite large range of 
the interaction strength if one assumes V{uj, Q) = V. If Ep ^ lOA, V ~ Ep, and A ^ 30meV , the correction to the 
position of the singularity, iVq ~ 20meV, is order of lO^^meV. Thus Xo is enough to describe the charge susceptibility 
as far as the most singular part is concerned. 

In the lattice, the Umklapp scattering is present and it leads to an additional contribution to the susceptibilty. This 
contribution provides another collective modes with the energy, ujq = 2E^q^ + 2/S.q^. As a result, the imaginary part 
of the susceptibility becomes 
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One can see that the shape of the structure factor depends on the ratio of the Fermi energy, Ep, and the maximum 
amplitude of the gap, A. In Fig. 2, we show the imaginary part of the charge susceptibihty for E'f/A = 7 and 2. It 
is usually assumed that the ratio of the Fermi energy and the gap is order of 5 ~ 10 because of the short coherence 
length though it has not been confirmed. In a recent tunneling experiment, it was claimed that the ratio might be 
order of one in the case of Bi2Sr2CaCu20s+y{BiSC'CO) [Q. As Ep/ A gets bigger, the anisotropy of the structure 
factor becomes stronger as shown in Fig. 2. If one could detect the charge collective modes by inelastic x-ray or 
electron scattering, it would provide the ratio of the Fermi energy and the gap of the cuprates. 
The charge susceptibility in the real space is found from the Fourier transform of Eq. (|^) as 

ImxC(c^, X, y) = —^^-—tlf!l[^^ry_j, (^^) + (2^ _ rl)J, {cor)], (12) 

where J,/(x) is the Bessel function and = + r^. Here r+ = {x + y) / {2^/2Ep) and r_ = (a; — ?/)/(2\/2A). Since 
the periodicity of the oscillation depends on the ratio Ep/A, the shape of the structure factor in the real space is 
enlongated along the r_ direction for Ep > A shown in Fig. 3. 

Let us now study the spin susceptibility for the momentum transfer near Q^r = (ti", tt). We found 
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Imx^(^, Q) = " ' 

It can be clearly seen in Eqs. and ( p^ that the spin susceptibility is smooth while the charge susceptibility diverges 
as mentioned above. The difference between the charge and spin susceptibilities comes from the coherence factors. 
The spin susceptibility has the maximum at Q,r = 0) for a given frequency lo and decreases as LOq approaches 
UJ. One can show that when ujq becomes larger than uj the imaginary part of the irreducible spin susceptibility is 

negligible. On the other hand, the real part of the irreducible spin susceptibility is {16 AEp)^-^ ^uj-^ — uj'' for ujq > uj. 

As a result, the spin susceptibility in the random phase approximation, x'^{^^ Q) = Xoi^j Q)/[l + Q)Xo('^i 
can support the collective modes if the interaction J(w,Q) — J > 16Ep. The shape of the spin structure factor 
depends on the ratio Ep/A. As mentioned above, the Umklapp scattering provides the additional contribution to 
the spin susceptibility. After including the Umklapp scattering, we find the following results. The positions of the 
magnetic peaks are near (tt, nil ± 5)) and (7r(l ± 5), tt), where 5 (in units of 27r/a) is given by 



1 + {16AEp/JY 

If < {IQAEp / J) and Ep/A 5 — 10, the peak splitting parameter, 5, is order of 0.1 and alomost independent of 
the frequency. In the experiments, 5 seems to be frequency independent Our results imply that the interaction 

strength, J, and the ratio of the Fermi energy and the amplitude of the gap determine the positions of the peaks for low 
frequencies. The weights of the magnetic peaks are also determined by the interaction strength. For Ep/A 1, the 
magnetic peaks form a circle with the radius S given by Eq. ( p^ ) which may happen in BiSCCO. However, it is known 
that the interaction is typically smaller than the Fermi energy in the superconducting state, so that the condition 
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for the interaction strength to support the collective mode is not realistic. Therefore, it suggests that either the 
charge fluctuations induce the anomalous interaction or another ingredient is needed to explain the incommensurate 
magnetic peaks found in the experiments. 

In the low energy limit, the susceptibilities for the momentum transfer near Qo = (0,0) are also interesting. It is 
found that 



Imxg(w,Q) 



1 oj^ - (V2Ag_) 
1 



Imxg(^, Q) - ^ : , (15) 
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where q = Q — Qo- Note that both the irreducible charge and spin susceptibilities are singular as lu approaches 
ujq. However, these singularities in the long wavelength limit are likely to be weakened by the vertex corrections. At 
realtively high frequencies ~ 2 A the charge and spin collective modes were previously found near Q = (2fci?,0) and 
{kp,kp) respectively. Q 

Using the above information, we compute the lowest order contribution to the NMR relaxation rate 1/Ti assuming 
that the scattering near dominates: 

1 [ d^Q Imxg(c.,Q) 

t;t = 1'^tJ^ — ^ — "(A^- 

Thus the relaxation rate l/Ti goes as which is consistent with the numerical calculations Q and an experiment 
for a range of the temperatures . 

We can also compute the quasiparticle scattering rate due to the charge fluctuations. We consider only the charge 
fluctuations because these are more stronger than those of the spin. The lowest order self-energy can be obtained 
from 

y = 1/2 f ^^'^ [ '^^ ^ Ck+QTs + Ak+QTi r dx Imxg(x,Q) 

J (2^)2 7 2^ + _ ^2^^ + J ^ u-x + ii^ ' ^ 

The imaginary part of the self energy in the / space is found at the node, ko, as 

POO p 

ImS/(a;, ko) Y J (2^^""^°^'^' ^^^"^ + + -^ko+q) + 6(0; - x - E^„+q)] 

-^"(^^ m 

It can be shown that the self energies in the ti and spaces have larger power in uj, O(aj^), compared to lu^. Thus, 
the leading contribution to the quasiparticle scattering rate, l/r(a;), is given by 

In summary, we presented the analytic results for the dynamic charge and spin susceptibilities near Qjr — (tt, tt) 
and Qo = (0, 0) in the BCS d-wave superconducting state. We found that the presence of the d-wave node leads to 
the incommensurate charge peaks for a set of wavevectors forming an ellipse near Qjr and Qo. We also showed that 
incommensurate magnetic peaks near Qtt can be obtained through the random phase approximation if the extremely 
strong interaction is assumed. When the Umklapp process is included, higher intensity of the magnetic structure 
factor would appear at (tt, 7r(l ± (5)) and (7r(l ± S),Tr) where the incommensurate magnetic peaks were found in the 
experiments on LSCO and YBCO [0j|J^. We also compute the NMR relaxation rate going as and the quasiparticle 
scattering rate as uj^ in the superconducting state. 

We thank G. Aeppli, S.-W. Cheong, P. Dai, H. Y. Hwang, C. M. Varma, and R. Walstedt for helpful discussions. We 
are especially greatful to Q. Si for explaining his previous and recent numerical results. One of authors (H.Y.K.) thank 
B. L. Altshuler and C. M. Varma for motivating the present problem and D. A. Bonn for explaining his experiment. 
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FIG. 1. The circular Fermi surface with a radius kF illustrating the wave verctor Q^. The diamond represents the tight 
binding Fermi surface at the half-filling. The node positions (ifc^/v^, ±fcF/v'2) in the case of the circular Fermi surface 
correspond to ko = (±7r/2, ±7r/2) in the tight binding model. As a result, the momentum transfer {■\/2kF, y/2kF) corresponds 
to = (7r,7r). 



FIG. 2. The charge susceptibility for (a) Ef/A — 7 and (b) 2 with the frequency 0.6 A. q+ and g_ are in units of 27r/a. 



FIG. 3. The charge susceptibility in the real space for Ef/A = 7 and the frequency 0.6A. x and y are in units of a. 
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